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ABSTRACT 

We re-examine MSW solutions of the solar neutrino problem in a two flavor scenario 
taking (a) the results on total rates and the electron energy spectrum from the 1117-day 
SuperKamiokande (SK) data and (b) those on total rates from the Chlorine and Gallium 
experiments. We find that the SMA solution gives the best fit to the total rates data from 
the different experiments. One new feature of our analysis is the use of the moments of the 
SK electron spectrum in a \ 2 analysis. The best-fit to the moments is broadly in agreement 
with that obtained from a direct fit to the spectrum data and prefers a Am 2 comparable to 
the SMA fit to the rates but the required mixing angle is larger. In the combined rate and 
spectrum analysis, apart from varying the normalization of the 8 B flux as a free parameter 
and determining its best-fit value we also obtain the best-fit parameters when correlations 
between the rates and the spectrum data are included and the normalization of the 8 B flux 
held fixed at its SSM value. We observe that the correlations between the rates and spectrum 
data are important and the goodness of fit worsens when these are included. In either case, 
the best-fit lies in the LMA region. 
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1 Introduction 



There is now strong evidence in support of an oscillatory behavior of neutrinos. The results 
on atmospheric neutrinos from SuperKamiokande |IJ find a comprehensive explanation in 
terms of oscillations, indicating a non-zero neutrino mass. This result and the solar neutrino 
problem, which has been held as a signal for neutrino oscillations for long, have created much 
excitement in the community. 

In this paper we examine the latest solar neutrino data from SuperKamiokande [|| along 
with those from the radiochemical Chlorine and Gallium experiments assuming that MSW 
resonant flavor conversion || is operative. The above experiments have presented the mea- 
sured arrival rates Q which, in all the cases, are less than the predictions from the Standard 
Solar Model (SSM). In addition, SuperKamiokande (SK) has provided the observed electron 
energy distribution [Q, ^] • We use these data to test the consistency of the MSW mechanism 
taken together with the SSM predictions in a two flavor scenario. There are several recent 
MSW analyses of the SK solar neutrino data || - [11]. Although the data fitting method fol- 
lowed in all of them is to minimize a x 2 function, the details of the statistical procedure used 
vary among the different groups. In this work, we indicate two different ways of performing 
the statistical analysis for the spectrum and the combined rate and spectrum data and check 
the consistency of the best-fit values of the mass squared differences and the mixing angles 
so obtained. For the analysis of the spectrum results we explore the possibility of using its 
moments as variables for fitting the data. The main advantage of these moments is that 
they probe the shape of the spectrum in a manner independent of the 8 B flux normaliza- 
tion uncertainties. For the combined analysis of the rates and the spectrum, apart from the 
standard procedure of varying the 8 B flux normalization and treating these two sets of data 
as independent, we also adopt a second method which takes into account the correlations 
among the rates and the spectrum data. 

In addition to the best-fit values of the oscillation parameters, we present the 99% C.L. 
and 90% C.L. allowed regions and the goodness of fit (g.o.f) of a particular solution. By g.o.f. 
is meant the probability that the x 2 will exceed Xmin- When presenting the allowed region 
we take Xmin t° be the value at the global minimum in that regionF]. For the neutrino fluxes 
and the neutrino production positions within the sun we use the BP98 solar model [12]. We 
consider oscillation of v e to a sequential (f„ or v T ) neutrino. 

This paper is organized as follows. In the next section we present the formulae for 
oscillation of neutrinos with the inclusion of matter effects both in the sun and during their 
passage through the earth. In section 3 we use the data on the total solar neutrino rates 
as measured at the Chlorine, Gallium, and SuperKamiokande (1117-day data) detectors to 
obtain the best-fit values of the neutrino mass splitting and the mixing angles. In section 4 
we consider the electron energy spectrum observed at SuperKamiokande. Using the MSW 
predictions, we obtain the best-fit values from a direct fit to the data as well as from a 
fit to the normalized moments. In section 5 we use both the total rates data and the SK 
electron energy spectrum data to make a combined fit. As noted earlier, here we allow 
the normalization of the 8 B spectrum to vary and compare the results with those obtained 



lr The other approach is to present the allowed regions with respect to the local minimum. 
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when the SSM prediction for this normalization is used, allowing the inclusion of correlations 
between the total rates and the observed spectrum via astrophysical uncertainties. We end 
in section 7 with a summary, discussions, and conclusions. 



2 Oscillation Probability 

In this work we restrict ourselves to the simplest case of mixing between two neutrino flavors. 
Assuming the neutrino mass eigenstates, Ui, reaching the earth to be incoherent [13|, the 
survival probabaility of an electron neutrino can be written as 

Pee = Pei^ie > (1) 
i 

where Pj is the probability of an electron neutrino state to transform into the i-th. mass 
state at the solar surface and is the conversion probability of the i-th mass state to the v e 
state after traversing the earth. One can express P ee in terms of the day-time (i.e. no earth 
matter effect) probability P® e asQ 

n (2P£-l)(sin 2 fl-P.g) 
^ e "^ ee+ cos 26 ' [) 

where 

p° = 0.5 + [0.5 - e(E u - E A )X]cos 26 M cos 26, (3) 

being the Heaviside function. 6 is the mixing angle in vacuum and 6m is the mixing angle 
in matter given by 

Am 2 sin 26 

tan 26 m = 5 7= • (4) 

Am 2 cos20- 2V2G F n e E u 

Here n e is the ambient electron density, E v the neutrino energy, and Am 2 (= m| — m 2 ) the 
mass squared difference in vacuum. 

E A = Am 2 cos 26/2V2G F n e \ pr , (5) 

gives the minimum v e energy that can encounter a resonance inside the sun, n e \ pr being the 
electron density at the point of production. X is the jump-probability between the mass 
eigenstates and for an exponential density profile, as is approximately the case in the sun, it 
is given by (Tf| 

x _ exp [-7T7fl(l - cos 26*)] - exp [-27rj R ] 
1 - exp[-27r7#] 

where 7r = 7 cos 26 / sin 2 26 and 

7r Am 2 sin 2 26 1 
4 E v cos26> 



' — a tt< „„„ no 1 d lnn e I ^ ' 



dr \ res 



2 This expression is not applicable for maximal mixing (cos 29 —0) [J14j and for this case we use eq. 
directly. 
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For calculating P^ e in eq. (|2|) we treat the earth as a slab of constant density (4.5 gm/cc). 
We have verified that this is a reasonable approximation since the location of the SK detector 
ensures that only in a rather small fraction of the time does the neutrino pass through the 
denser core. 

The definition of x 2 used by us in the following sections is, 

x 2 = e { p t - p r) (^ 2 ) (if - F e r) ■ w 

Here Ff = T$ /T^ P98 where £ is th (for the theoretical prediction with oscillations) or exp (for 
the experimental value) and T{ stand for the quantities being fit (total rates from different 
experiments, electron energy spectrum for different energy bins, etc.). The error matrix cr™ 
contains the experimental errors, the theoretical errors and their correlations. 

3 Total rates 

In this section we perform an analysis of the total rates as measured at the various experi- 
ments. The data that we use for the total rates are given in Table 1. In particular, we use the 
1117-day SK data. For the Ga experiments we take the weighted average of the SAGE and 
Gallex results. Because SK has better statistics we do not include the Kamiokande results. 



Experiment 


Chlorine 


Gallium 


SuperKamiokande 


Observed Hate 
BP98 Prediction 


0.33 ± 0.029 


0.562 ± 0.043 


0.465 ±0.015 



Table 1: The ratio of the observed solar neutrino rates to the corresponding BP98 SSM 
predictions used in this analysis. The results are from Refs. [f|] and [gj. For Gallium, 
the weighted average of the SAGE and Gallex results has been used. 



To get the best-fit values of the parameters, we minimise a x 2 function defined as in eq. 
(||). For evaluating the error matrix, <jy, we use the procedure described in [jn|. Our best-fit 
results for the total rates are summarized in Table 2. These fits have 1 degree of freedom (3 
experimental data points - 2 parameters). It is clear from this Table that the small mixing 
angle solution for the sequential neutrino case gives by far the best fit to the total rates data. 

It is seen from Table 2 that the best-fit values obtained for the 1117-day data are not 
markedly different from those obtained for the 825-day data. This indicates that the best-fit 
points obtained from the analysis of total rates are quite robust and are not expected to 
change drastically with more data from SK. However, the quality of fit is quite sensitive to 
these small changes: the g.o.f. has become a little poorer for the SMA solution while for both 
the LMA and LOW solutions it has improved with the accumulation of more SK data. 

In Fig. 1 we show the 99% C.L. and 90% C.L. contours in the sin 2 9 — Am? plane for the 
MSW solution. There is no solution in the dark side (sin 2 9 > 0.5) from the analysis of total 
rates. Note that in [17| allowed regions were obtained in the dark side from the total rates 
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No. of days of 
SK running 


Nature of 
Solution 


Am 2 
in eV 2 


sin 2 9 


/V. 772 % 71 


Goodness 
of fit 


1117 


SMA 


5.19 x 10"° 


1.53 x KT 3 


0.30 


58.38% 


LMA 


1.42 x 1(T 5 


0.22 


3.52 


6.06% 


LOW 


9.97 x 10-* 


0.38 


6.64 


.997% 


825 


SMA 


5.33 x 10 _e 


1.50 x 10~ 3 


0.06 


80.64% 


LMA 


1.40 x 10"° 


0.23 


4.31 


3.78% 


LOW 


9.98 x 1(T 8 


0.38 


7.39 


0.65% 



Table 2: The best-fit values of the parameters, Xmim an d the g.o.f. for fits to the total rates 
of the different experiments. SMA, LMA, LOW stand for the Small Mixing Angle, 
Large Mixing Angle and Low mass-low x 2 regions, respectively. The results for the 
earlier 825-day SK data (0.475 ± 0.015) are also presented for comparison. 



analysis while in [18| no allowed region was found for > ir/4. We agree with [|J|; a possible 
origin of this is that the same numerical density profile of the sun from BP98 has been used 
in these analyses. Our results for the total rates for the 1117 day data are in agreement with 
the analysis given in [19]. 

4 Observed spectrum and its moments 

In addition to the total rates, SK has provided the number of events (normalized to the SSM 
prediction) in 18 electron recoil energy bins of width 0.5 MeV in the range 5.0 MeV to 14 
MeV and a 19th bin which covers the events in the range 14 to 20 MeV Q] . The systematics 
of the first bin are still under study and for our analysis we do not use it. 

4.1 Observed spectrum 

In this subsection we present the results obtained by directly fitting the SK spectral data. 
The theoretical predictions are calculated bin by bin and in the fitting procedure, in addition 
to the neutrino mixing parameters Am 2 and sin 2 9, we also allow the absolute normalization 
of the 8 B flux, Xb, to vary[]. The error matrix aij used by us (see eq. (||)) is Q 

(?ij)sp = fiiji&i^tat ®i,uncorr) + ®i,expO "j,exp + ^ifial^j^cali (9) 

where we have included the statistical error, the uncorrelated systematic errors and the 



energy-bin-correlated experimental errors [20] as well as those from the calculation of the 
shape of the expected spectrum pi]] . Since we vary the normalization of the 8 B flux we do 
not include its astrophysical uncertainties separately. 
The best-fit point from this analysis is found to be 



J X B = 1, for the SSM. 
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• Am 2 = 2.29 x 1(T 6 eV 2 , sin 2 9 = 0.009, X B = 1.4, x min = 9 - 46 > g-°- f - = 85.23%. 

For these values of Am 2 and sin 2 9, choosing Xb = 1, the data on total rates give a bad 
fit ( x 2 = 73.82) as the high energy 8 B neutrinos get suppressed more than observed. For the 
spectrum fit this problem can be avoided by a high Xb > 1. Since the 8 B flux normalization 
is allowed to vary, a large range of Am 2 and sin 2 9 including the 9 > it/A region remains 
allowed by the spectral data as is shown in Fig. 2. Only the area inside the contours is 
disallowed at 90% C.L. from spectral data analysis. 

To further examine this mismatch between the fits to the total rates and those to the SK 
spectral data, in Fig. 3 we present the behavior of x 2 as obtained from the spectral analysis 
if we keep the parameters sin 9 and Am 2 in the SMA, LMA, and LOW regions of the fit 
to the total rates. In Figs. 3(a), 3(b), and 3(c) are shown the variation of the x 2 from the 
spectral data with sin 2 9 lying in the SMA, LMA, and LOW regions respectively. For any 
chosen sin 2 9, we let Am 2 vary over the corresponding range permitted by the fit to the total 
rates at 99% C.L. (from Fig. 1) and plot the minimum value found. We consider two cases, 

• the 8 B normalization is held fixed at its SSM value (solid curves) 

• the 8 B normalization is permitted to vary (broken curves). 

For Figs. 3(d), 3(e), and 3(f), the roles of sin 2 9 and Am 2 are interchanged. Fig. 3 indicates 
that if we allow the 8 B normalization to vary then the Am 2 and sin 2 9 allowed at 99% C.L. 
from the total rates are allowed at 90% C.L. from the spectral analysis. 

Figs. 2 and 3 lead us to the conclusion that, at the present moment, the SK electron 
spectral data do not provide tight controls over the allowed parameter range. 



4.2 Moments of spectrum 

Since the absolute normalization of the 8 B flux is not precisely known, it is of interest to 
look for variables which probe neutrino oscillation effects in the data in a manner immune 
to this uncertainty. Normalized moments of the observed electron spectrum can be useful as 
one such set of variables |2^, |2^]. In practice, to compare with the data, it is convenient to 
standardize with respect to the SSM predictions by using 



Mr, 









_{N(Ei)\sSM _ 













(10) 



where E{ is the mean energy of the i-th bin and N{Ei) is the number of events in this bin. 
Depending on whether the experimental or the theoretically predicted value of the variable 
is under consideration, N(Ej) is obtained either from experiments or from the theoretical 
model under test. It is clear that these variables carry information about the shape of the 
neutrino spectrum which, if oscillations are operative, undergoes modification from the SSM 
prediction due to the energy dependence of the survival probability. It is obvious that the 
above moments are independent of the absolute normalization of the 8 B flux. 

We have calculated the moments of the 1117-day data on the electron energy spectrum 
presented by SuperKamiokande 0. These are presented in Table 3. The error in the higher 
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moments increases rapidly with the order and the ones beyond the sixth are not of much use. 



Order of 


Value of 


Calculated Error 


Moment 


Moment 


in Moment 


1 


10.25 


0.33 


2 


113.84 


7.12 


3 


1356.59 


164.34 


4 


17167.7 


3616.65 


5 


228876.2 


75981.0 


6 


3188064.0 


1541763.0 



Table 3: Moments of the observed electron energy spectrum and their calculated errors 
obtained from the SK (1117 days) data. 

Using these variables in a x 2 - anai ysis we And that fitting the first four moments results 
in the same best-fit values of Am 2 and sin 2 6 as those obtained using the first five or six 
moments. The best-fit region^ is found to be 

• Am 2 = 6.94 x lO" 6 eV 2 , sin 2 9 = 0.007 - 0.010, x 2 min = 0.0001, g.o.f = 99.995%. 

The very small value of Xmin f° r this fit should not be regarded as a major success of the 
theory but rather reflects the large errors associated with the moments as obtained from the 
present data. It is gratifying that the best-fit values obtained by this method are in the same 
broad region as those from fitting the recoil electron energy spectrum. 



5 Combined fits to rates and spectrum 

In this section we present the results of the combined fit to the total rates and the spectrum 
data. We have performed this global fit by the following two methods: 

(a) We treat the rates and the electron spectrum data as independent. In this approach 
we vary the 8 B flux normalization as a free parameter. 

(b) We fix the 8 B flux normalization at the SSM value (=1) and include the correlations of 
the 8 B flux uncertainty between the rates and spectrum data. To our knowledge, this 
approach has not been pursued in any previous analysis. 



4 The x 2 remains unchanged when sin 2 8 is varied over the indicated range 
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5.1 Fits using the B flux normalization as a free parameter 

For this case the definition of \ 2 is, 

x 2 = E ( F t - F i xp ) K?) (*f - F ! xp ) 

+ E (x B ^ h -Rr P )(^ 2 )s P {xBRf-RT P ), (11) 

ij=l,18 

where the first term on the r.h.s is from the fit to the total rates and the second from that to 
the spectral data. As we allow the normalization of the 8 B flux to vary as a free parameter 
we switch off the SSM astrophysical uncertainties arising because of this component. Since 
it is the 8 B flux that enters the rates as well as the spectrum data, in this manner of fitting 
the data the correlations between the rates and the spectrum are absent; the error matrix is 
block diagonal and one can treat Xrate an d xlpectrum as independent. There are 18 (= 21 - 
3) degrees of freedom in this case. The best-fit values we obtain are presented in Table 4. 



Nature of 
Solution 


Am 2 
in eV 2 


sin 2 


X B 


x 2 ■ 


Goodness 
of fit 


SMA 


4.81 x l(T e 


5.92 x 10~ 4 


0.61 


12.62 


81.36% 


LMA 


1.82 x 1(T 5 


0.18 


1.39 


10.97 


89.56% 


LOW 


1.0 x 1(T Y 


0.38 


0.95 


17.21 


50.87% 



Table 4: The best-fit values of the parameters, Xmiw an d the g.o.f. for fits to the total rates 
measured at the CI, Ga, and SK detectors along with the electron energy spectrum 
from SK (1117 days) when the 8 B flux normalization factor Xb is allowed to vary. 

In Fig. 4a we show the 99% and 90% C.L. allowed regions for the combined analysis of 
total rates and the observed electron spectrum. The best-fit points in these plots are obtained 
by varying Xb in addition to Am 2 and sin 2 0. 

5.2 Fits including correlations between rates and spectrum via 8 B flux 

For this case we include the correlations between the theory errors in the rate and the 
spectrum data. This comes through the 8 B flux, as it enters both. Since we include the 
SSM astrophysical uncertainties in the 8 B flux the normalization factor for it is held fixed at 
the SSM value. Now the individual x 2 due to the spectrum and the rates cannot be summed 
independently and the combined \ 2 is defined as, 

X 2 = E ( F i h - F i XP ) (<^ 2 ) ( F j k - F j XP ) ' (12) 

M=l,21 

where the G{j is now a 21 x 21 matrix defined in the following way, 
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For i,j = 1 ... 3 

a ij = ( a ij)th + (Oij)exp, (13) 

where 

8 8 ell 

(o-y)tfi = <*ij I] i& (AC ai ) 2 + ^ i? ai a *k apk (A In X k ) 2 . (14) 

a=l o,/3=l k=l 

where the first term is due to the cross-section uncertainties and the second term (c ap ) 
is due to the astrophysical uncertainties [16]. The off-diagonal elements in the error 
matrix come through a ap . R a i denotes the contribution of the a-th source to the rate 
of the i-th experiment. a ak = 51n(/) a /51nX k , where 5lncf) a is the error in the a-th 
component of the spectrum due to the input parameter X k . 

For i = 4 ... 21 and j = 1 ... 3 

8 11 

ofj = R8 BiRaj as Bk a ak (/^\uX k ) 2 . (15) 
a=l k=l 

For i = 1 ... 3 and j = 4 ... 21 



8 11 

afj = Y RaiRsBj Y a a ka8 Bk (AlnX k ) 2 . (16) 

a=l k=l 



for i = 4 ... 21 and j = 4 ... 21 

ii 



°fj = ( a ij)sp + R»Bi R »Bj Y as Bk as Bk(A l n ^fc) 2 • (17) 



k=l 

2 



In this case the number of degrees of freedom is 19 (= 21-2). The Xmin an d l ne best-fit 
values we obtain are shown in Table 5. 



Nature of 
Solution 


Am' 2 
in eV 2 


sin 2 9 


x 2 ■ 


Goodness 
of fit 


SMA 


5.15 x lO" 6 


4.8 x 10" 4 


15.72 


67.58% 


LMA 


2.23 x 1(T 5 


0.25 


14.53 


75.19% 


LOW 


1.00 x 10~ Y 


0.38 


17.62 


54.79% 



Table 5: The best-fit values of the parameters, Xmiw an( l * ne §-°-f- f° r to the total rates 
measured at the CI, Ga, and SK detectors along with the electron energy spectrum from 
SK (1117 days). In this case the 8 B flux is chosen as in the SSM and the correlation 
between the rates and the spectrum due to astrophysical uncertainties is included. 
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We find that the fits are of poorer quality in this case as compared to the previous one. 

In Fig. 4b we show the 99% and 90% C.L. allowed regions for the combined analysis 
of total rates and the observed electron spectrum for MSW conversion to sequential neutri- 
nos including the correlations between the rates and the spectrum due to the astrophysical 
uncertainties of the 8 B flux normalization. 



6 Summary, Discussions, and Conclusions 

In this paper we have performed a detailed x 2_ana lysis of the latest SK solar neutrino data 
together with the results from the CI and Ga experiments in terms of two-generation MSW 
conversions of v e to sequential (v^, u T ) neutrinos. 



Compared to the recent analyses in the literature ||-[11| there are two new features in 
our analysis. 

• We fit the observed electron energy spectrum data in two different ways, exploring for 
the first time, the use of moments of the energy spectrum in a x 2_ana lysis. 

• The combined fits to the total rates and spectral data are also performed in two different 
manners. In the first, the 8 B flux normalization is used as a free parameter while in 
the other the SSM normalization is chosen for it and correlations between the rates and 
spectral data due to astrophysical uncertainties of the 8 B flux are included. 

We find that the two- generation MSW scenario can well explain the data on total rates. 
The solution in the SMA (Small Mixing Angle) region is preferred over the other possibilities 
although the quality of the fit is poorer as compared to the one obtained using the 825-day 
SK data. 

The best-fit from the spectrum data comes in a region disallowed from the total rates. In 
this region the 8 B neutrinos are suppressed much more than required by the rates data. For 
the analysis of the spectrum, the absolute normalization of the 8 B flux, Xb, has been permit- 
ted to be greater than unity, thus effectively compensating the shortfall. We have explored 
the use of normalized moments of the observed electron energy spectrum to signal MSW 
resonant flavour conversion. These variables are independent of the absolute normalization 
of the 8 B flux and probe the effect of oscillations on the spectral shape. This procedure is 
somewhat handicapped by the large errors on the moments calculated from the present data. 
However, the best-fit values obtained by the two methods are more or less in agreement. 

Similarly, for the two methods followed in the combined x 2 analysis of the rates and time 
averaged spectrum data, the best-fit values are not much different. The first approach gives 
a better fit because we utilise the freedom of varying the 8 B flux normalization. We remark 
that in the combined analysis, where the 8 B normalization is held fixed at the SSM value, 
the correlations between the rates and the spectrum data are found to be important and 
thus one should use caution regarding results obtained treating these as independent. For 
both methods, the best-fit from the combined analysis falls in the LMA region. Compared 
to the rates analysis the goodness of fit of the LOW(SMA) region increases(decreases). With 
the inclusion of the day-night dependence of the data the goodness of fit in the SMA region 
worsens further 101. 
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In this work we have not included the new GNO result [E4j] which is consistent with 
the Gallex and SAGE data. Thus its inclusion is not expected to affect the conclusions 
drastically. For illustration we give below the results of the global analysis of rates and 
spectrum including the GNO data. We take the weighted average of Gallex and GNO and 
treat SAGE as a separate experiment. The best-fit values and Xmin ^ na * we are: 

• sin 2 9 = 5.26 x 10" 4 , Am 2 = 5.28 x 10" 6 eV 2 , X B = 0.61, x 2 min = 12 - 73 > g-°- f = 
85.21% (SMA) 

• sin 2 6> = 0.18, Am 2 = 2.48 x 10~ 5 eV 2 , X B = 1.39, xLm = 1L55 > §-°- f = 
90.39% (LMA) 

• sin 2 6> = 0.41, Am 2 = 9.39 x 10" 8 eV 2 , X B = 0.89, xLn = 19 - 85 > g-°- f = 
40.34% (LOW) 

Thus the global best-fit continues to be in the LMA region. 

In conclusion, we have probed the most recent solar neutrino data on total rates and 
the observed electron energy spectrum at SK from various angles within the framework of 
MSW flavour conversion. We find good fits in some situations but a degree of uncertainty 
still remains since different fits do not prefer the same values of the parameters. More data 
from the running and new experiments, it is hoped, will further sharpen the results in the 
near future. 
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Figure Captions 



Fig 1. The 99% and 90% C.L. allowed regions in the Am 2 - sin 2 9 plane from the analysis 
of total rates for the Chlorine and Gallium detectors and the 1117-day data from SK. The 
best-fit points are also indicated. The dark side (sin 2 9 > 0.5) is indicated by the dashed line. 

Fig 2. The 90% C.L. allowed region in the Am 2 - sin 2 9 plane from the 1117-day SK recoil 
electron spectrum data. The regions enclosed by the contours are disallowed. The best-fit 
point is indicated. The dark side is to the right of the dashed line. 

Fig. 3. The minimum x 2 for fits to the SK 1117-day recoil electron spectrum as a function of 
sin 2 ^ (Am 2 ) are shown in (a), (b), and (c) ((d), (e), and (f)) when the parameter ranges are 
determined by the 99% C.L. allowed regions in the SMA, LMA, and LOW fits respectively to 
the total rates data. The solid (broken) curves are obtained when Xb is held fixed at its SSM 
value (allowed to vary). The dash-dotted line indicates the 90% C.L. limit for 3 parameters. 
See text for more details. 

Fig. 4. The 99% and 90% C.L. allowed region in the Am 2 - sin 2 9 plane from an analysis 
of the total rates from the Chlorine and Gallium detectors and the 1117-day SK data taken 
together with the 1117-day SK recoil electron spectrum. The normalization of the 8 B flux 
is chosen as a free parameter in (a) and held fixed at the SSM value in (b). In (b) the 
correlations between the rates and spectrum data are included. The best-fit points are also 
indicated. The dark side corresponds to the right of the dashed line. 
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